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Atomic-layer  2D  crystals  have  unique  properties  that  can  be  significantly  modified  through  interaction 
with  an  underlying  support.  For  epitaxial  graphene  on  SiC(OOOl),  the  interface  strongly  influences  the 
electronic  properties  of  the  overlaying  graphene.  We  demonstrate  a  novel  combination  of  x-ray  scattering 
and  spectroscopy  for  studying  the  complexities  of  such  a  buried  interface  structure.  This  approach 
employs  x-ray  standing  wave-excited  photoelectron  spectroscopy  in  conjunction  with  x-ray  reflectivity  to 
produce  a  highly  resolved  chemically  sensitive  atomic  profile  for  the  terminal  substrate  bilayers,  interface, 


and  graphene  layers  along  the  SiC[0001]  direction. 
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Epitaxial  graphene  (EG)  grown  on  the  Si-terminated 
face  of  silicon  carbide  [SiC(OOOl)]  is  of  interest  due  to 
its  unique  physical  properties  [1,2],  potential  in  numerous 
applications  [3-5],  and  amenity  for  wafer-scale  fabrication 
[6,7].  EG/SiC(0001)  production  proceeds  via  preferential 
Si  sublimation  and  reconstruction  of  the  surface  C  to  yield 
a  heterostructure  that  exhibits  the  technologically  relevant 
properties  characteristic  of  graphene,  such  as  the  Dirac 
cone  band  structure  [8,9]  and  large  room-temperature 
Hall  mobilities  [1,6].  As  such,  EG/SiC(0001)  is  a  com¬ 
paratively  viable  option  for  the  large-scale  production  of 
graphene.  However,  despite  the  advancement  in  both  fun¬ 
damental  understanding  and  practical  application  of  the 
material,  there  remains  much  unknown  about  the  structure 
and  behavior  of  EG/SiC(0001),  particularly  concerning  the 
influence  of  the  EG/SiC(0001)  interface  on  the  properties 
of  the  overlaying  graphene. 

Early  studies  revealed  that  EG/SiC(0001)  possesses  a 
complex  6^3  X  6^3R30°  (6R3)  reconstructed  interfacial 
layer  [10],  referred  to  herein  as  the  interfacial,  or  EGq, 
layer.  This  layer  has  significant  influence  on  the  growth, 
morphology,  and  electronic  behavior  of  the  overlaying 
graphene  [8,11-14].  It  has  been  contentiously  debated  as 
the  cause  of  EG/SiC(0001)  symmetry  breaking  and  band 
gap  opening  [15,16],  and  has  been  identified  as  a  primary 
source  for  carrier  doping  and  scattering  in  the  graphene 
layers  [17,18].  Work  on  H-intercalated  graphene  transis¬ 
tors  highlights  the  influence  of  this  layer  on  the  electronic 
properties  of  the  overlaying  graphene,  demonstrating  that 
physically  decoupling  the  EGq  layer  from  the  substrate 
improves  properties  critical  to  graphene-based  transistor 
performance  [19,20].  A  precise  understanding  of  the  struc¬ 
ture  of  EG/SiC(0001)  is  therefore  not  only  fundamentally 
important,  but  is  also  central  to  the  prospects  of  graphene 
engineering  and  functionalization. 


PACS  numbers:  61.48.Gh,  61.05.cm,  68.49.Uv,  79.60.-i 

Because  of  the  importance  of  the  interfacial  layer  to  the 
behavior  of  EG/SiC(0001),  there  have  been  numerous  efforts 
to  characterize  its  structure,  including  low-energy  electron 
diffraction  [10,13],  scanning-tunneling  microscopy,  [21-23], 
x-ray  photoemission  spectroscopy  (XPS)  [24,25],  x-ray 
reflectivity  (XRR)  [26],  and  cross-sectional  transmission 
electron  microscopy  [27,28].  These  experimental  reports, 
combined  with  computational  studies  [14,29,30],  have  gen¬ 
erated  a  multitude  of  interfacial  models.  These  include  those 
with  weak  [10]  or  strong  [12,14,25,29]  substrate-interlayer 
coupling,  C-rich  [14,23,25]  or  Si-rich  interfacial  structures 
[21,26],  or  significant  populations  of  Si  dangling  bonds 
[12,29].  However,  despite  this  persistent  incongruence,  it  is 
now  understood  that  the  interface  consists  of  a  corrugated, 
topologically  graphenehke  layer  which,  due  to  its  coupling 
with  the  underlying  SiC,  lacks  the  electronic  properties  of 
pristine  graphene  [25,31]. 

In  this  Letter  we  detail  the  structure  of  the  interface  by 
employing  a  suite  of  x-ray  characterization  techniques, 
including  depth- sensitive  XPS,  x-ray  standing  wave- 
enhanced  XPS  (XSW-XPS),  and  x-ray  reflectivity  (XRR). 
These  tools,  when  employed  collectively,  provide  the 
chemically  specific  structural  information  necessary  to 
clarify  previously  unknown  details  of  the  EG/SiC(0001) 
interface.  This  approach  ultimately  enables  the  construc¬ 
tion  of  a  chemically  resolved  interfacial  map  with  sub-A 
resolution  along  the  SiC [0001]  direction. 

The  XSW  technique  affords  conventional  photoelectron 
spectroscopy  with  high  spatial  resolution  due  to  the  influ¬ 
ence  of  the  XSW  [here  produced  by  the  SiC(0006)  Bragg 
reflection]  on  the  photoabsorption  process.  A  depiction  of 
this  phenomenon  is  shown  in  Pig.  1 .  The  XSW  phase  can  be 
adjusted  relative  to  the  atomic  planes  by  tuning  the  incident 
beam  energy  (E^)  within  the  range  of  the  Bragg  reflection, 
which  in  turn  modulates  the  x-ray  photoelectron  yields  of 
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FIG.  1  (color  online).  Schematic  of  the  XSW  generated  by  the 
6//-SiC(0006)  Bragg  reflection  (period  of  d  =  2.52  A).  The 
incident  and  reflected  plane  waves  (wave  vectors  Rq  and  kn, 
respectively)  interfere  to  produce  a  standing  wavefield  (antin¬ 
odes  shaded  red). 


atoms  within  the  field.  Depending  on  the  spatial  distribution 
of  these  atoms  with  respect  to  the  substrate,  this  process 
produces  distinct  photoemission  yield  modulations  as  a 
function  of  Ey.  XSW  data  analysis  determines  two  model- 
independent  fit  parameters,  fj  and  Py,  which  are  the  ampli¬ 
tude  and  phase  of  the  H  =  0006  Eourier  component  of  the 
jth  chemically  sensitive  atomic  density  profile.  Eor  details, 
see  the  Supplemental  Material  [32]  and  Ref.  [33]. 

In  this  work  we  will  assume  a  normalized  Gaussian 
distribution  for  each  XPS-selected  atomic  layer.  The 
(0006)  Eourier  component  is  therefore 

Fj  = 

Eor  this  model,  the  XSW-measured  Eourier  amplitude  fj 
determines  the  Gaussian  width  cr^,  and  the  measured  phase, 
Pj  =  Zj Id  (0  <  Py  <  1),  determines  the  mean  height  Zj 
(modulo-d)  of  the  jth  atomic  layer  relative  to  the  substrate 
unit  cell  origin  zq  (Eig.  1).  Because  of  the  modulo-d  ambi¬ 
guity  and  potential  contributions  to  Pj  from  multiple  atomic 
layers,  we  supplement  this  analysis  with  high-resolution 
XRR  (see  Ref.  [34]  and  references  therein),  a  measurement 
that  is  sensitive  to  electron  density  distributions  over  a 
larger  spatial  extent. 

EG/SiC(0001)  samples  of  1.3  and  0.5  monolayer  (ML  of 
graphene  =  38.2  C/nm^)  graphene  coverage  (@7^)  were 
prepared  by  thermal  decomposition  of  6//-SiC(0001)  sub¬ 
strates  in  ultrahigh  vacuum  (UHV)  at  1200-1300  °C. 
Results  for  the  1.3  ML  sample  are  presented  in  the  main 
text.  Complimentary  analysis  of  the  0.5  ML  sample,  as  well 
as  an  Ar-grown  sample  of  1 .7  ML  coverage,  are  provided  in 
the  Supplemental  Material  [32]  to  confirm  and  generalize 
the  conclusions  of  the  main  text.  Prior  to  measurements,  the 
samples  were  annealed  in  UHV  at  —600  °C,  and  the  pres¬ 
ence  of  the  6P3  interface  was  verified  with  low-energy 
electron  diffraction.  Eurther  sample  preparation  details  are 
provided  in  the  Supplemental  Material  [32]. 

XSW-XPS  measurements  were  performed  in  UHV  at 
the  ID32  beam  line  of  the  European  Synchrotron 
Radiation  Eacility.  Working  near  the  SiC(0006)  Bragg 


back-reflection  condition  {6^  —  88°),  the  incident  beam 
energy  and  EWHM  bandwidth  were  Ey  =  =  2450 

and  ^.Ey  =  0.34  eV,  respectively.  The  total  EWHM  reso¬ 
lution  of  the  spectrometer  was  0.6  eV.  To  vary  the  depth 
sensitivity,  the  photoemission  angle  {a)  was  set  to  either 
78°  or  2°,  as  shown  in  Eig.  2.  The  grazing  (of  =  2°) 
emission  geometry  was  used  during  XSW  experiments. 
XRR  data  were  acquired  at  5ID-C  at  the  Advanced 
Photon  Source.  XRR  integrated  signal  intensity  was 
extracted  following  the  methods  described  in  Ref.  [34]. 
XRR  data  fitting  was  limited  to  the  region  1.2  <  L  <  14.0, 
where  L  =  / {lir)  is  the  SiC  reciprocal  lattice 

index,  =  Air  sin(2^/2)/A  is  the  out-of-plane  component 
of  the  momentum  transfer  vector,  26  is  the  scattering 
angle,  A  is  the  x-ray  wavelength,  and  c^ic  =  15.12  A  is 
the  lattice  constant  for  6//-SiC. 

Eigure  2  presents  depth- sensitive  (effective  probing 
depth  Agff)  core-level  photoelectron  spectra  used  to  iden¬ 
tify  near-surface  C  and  Si  chemical  species.  When  mea¬ 
sured  at  emission  angle  a  =  1S°  [Eigs.  2(a)  and  2(b)],  the 
Si  Iv  and  C  Is  spectra  are  dominated  by  bulk  SiC  compo¬ 
nents  with  binding  energies  (BE)  of  1841.6  ±0.1  and 
283.8  ±0.1  eV,  respectively.  The  C  Is  spectrum  possesses 
additional  high-energy  core-shifted  components  from  gra¬ 
phene  and  interfacial  layers,  as  previously  investigated 
[22,24,25].  The  Si  Iv  spectrum,  in  contrast,  appears  to 
consist  of  a  single  component.  By  reducing  a  to  2° ,  A^ff 
is  reduced  from  nm  to  A  scale  [Eigs.  2(c)  and  2(d)]. 
A  comparison  of  the  C  Is  spectra  in  Eigs.  2(a)  and  2(c) 
reveals  strong  attenuation  (X  1/30)  of  the  Csuik  integrated 
signal  with  respect  to  the  higher  BE  spectral  components. 


Binding  Energy  (eV) 


1845.0  1842.5  1840.0 


Si  1s  ^ 


FIG.  2  (color  online).  C  Is  and  Si  U  photoelectron  spectra  for 
1.3  ML  EG/SiC(0001).  Sisuik^  Csuik^  EG,  Si,  and  S2  components 
are  shown  in  purple,  blue,  green,  red,  and  brown,  respectively. 
(a)-(b)  Data  collected  at  emission  angle  a  =  78°.  (c)-(d)  Data 
collected  with  a  highly  surface-sensitive  geometry  (a  =  2°). 
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In  contrast,  no  distinguishable  surface- specific  component 
is  present  in  the  of  =  2°  Si  spectrum  [Fig.  2(d)]. 

We  identify  three  distinct  C  surface  species,  labeled 
in  Figs.  2(a)  and  2(c)  as  EG,  Si,  and  82-  These  components 
are  shifted  in  BE  by  +0.95,  +1.35,  and  +1.95  eV  with 
respect  to  the  Csuik  peak,  and  have  resolution-broadened 
EWHMs  of  0.7,  1.0,  and  1.1  eV.  The  EG  component  arises 
from  the  graphene  layers,  and  the  Si  and  S2  components 
belong  to  species  within  the  interfacial  layer.  The  Si  1^ 
spectra  are  fit  with  a  single,  slightly  asymmetrically  broad¬ 
ened  line  shape  [Eigs.  2(b)  and  2(d)].  When  the  geometry 
is  more  highly  surface  sensitive  the  peak  broadens  by 
—20%,  as  would  be  consistent  with  the  increased  relative 
contribution  of  strained,  near-surface  Si-C  bonds  possess¬ 
ing  a  variety  of  bonding  configurations  due  to  interactions 
between  terminal  Si  species  and  C  within  the  interfacial 
layer. 

Using  the  XPS  peak-fitting  models  described  above, 
we  extract  the  chemically  specific,  E^-dependent  XSW 
photoemission  yields  (Eig.  3).  Each  of  the  five  measured 
components  exhibits  a  distinct  XSW  modulation,  demon¬ 
strating  that  each  chemical  species  possesses  a  distinct 
distribution  profile  along  the  SiC[0001]  direction.  Results 
from  least-squares  fitting  of  the  yield  equation  [32,33]  to 
the  data  are  overlaid  for  each  component  in  Eig.  3.  All  best- 
fit  fj  and  Pj  values  are  reported  in  Table  I. 

The  Eourier  phases  for  both  Sisuik  ^Buik  indicate  that 
these  species  are  centered  at  their  nominal  SiC  positions 
(equal  to  0.00  and  —0.25,  respectively).  However,  small 
deviations  of  XSW  best-fit  results  for  both  species  from  the 
ideal  Sisuik  nnd  values  (—0.05  and  +0.05  A,  respec¬ 
tively)  imply  the  possibility  of  small  displacements.  The 
Eourier  amplitude  =  1.0  ±  0.1  denotes  essentially 
perfect  Csuik  coherency  with  the  substrate.  Conversely, 
the  SiBuik  result  =  0.8  ±  0.1)  suggests  a  broader- 

than-ideal  distribution  for  Si  atoms  within  the  topmost  SiC 


Ey-’Es  E^-Es 


FIG.  3  (color  online).  XSW  results  for  (a)  near-surface  bulk 
and  (b)  surface  species  corresponding  to  XPS  components  in 
Figs.  2(c)  and  2(d).  Measured  (gray  circles)  and  fitted  (black) 
reflectivity  from  the  SiC(0006)  Bragg  peak  are  plotted  vs  energy 
offset  from  the  geometrical  Bragg  condition  {E^  =  2450  eV). 
Normalized  yields  and  best-fit  results  (black)  are  offset  for  clarity, 
and  simulated  yields  for  ideal  SiC  bulk  are  in  dashed  black. 


bilayers  =  0.3  ±  0.1  A,  as  compared  to  the  bulk 

—0.1  A  thermal  vibrational  amplitude  for  SiC  at  room 
temperature). 

The  results  for  the  Si  and  S2  components  are  shown  in 
Fig.  3(b).  Converted  to  absolute  positions  [Eq.  (1)],  the  Si 
and  S2  components  reside  at  zs^  =  2.39  ±  0.13  and  zs^  ^ 
2.07  ±  0.10  A  above  the  terminal  Si  layer.  In  addition,  a 
comparison  of  the  Eourier  amplitudes  of  the  interfacial 
species  reveals  a  smaller  value  for  Si,  indicating  that 
is  significantly  broader  in  comparison  to 

The  remaining  C  Is  component  arises  from  atoms 
located  in  the  EG  overlayers.  Because  EG/SiC(0001) 
does  not  grow  in  a  single,  uniform  monolayer  [6],  single¬ 
layer  XSW  modeling  is  inappropriate.  Instead,  the  /eg  and 
Peg  values  are  determined  from  the  superposition  of  con¬ 
tributions  from  k>  I  atomic  layers  of  graphene.  Here,  the 
long-range  structural  information  lacking  in  XSW  is  pro¬ 
vided  by  XRR,  thereby  enabling  the  exploration  of  models 
with  varying  EG  layer  coverage. 

XRR  data  are  analyzed  by  comparing  measured  reflec¬ 
tivity  values  to  those  calculated  from  a  model  consisting  of 
k  atomic  layers,  each  described  by  layer  occupancy  (S^), 
position  (zk),  and  distribution  width  (cr^)  [26,34].  In  this 
work,  we  are  able  to  constrain  the  XRR  analysis  using 
XSW  results  (within  1 -sigma  confidence  levels),  thereby 
mitigating  ambiguities  that  commonly  arise  during  model- 
based  XRR  fitting.  Eor  SiBuik  and  CBuik^  we  allow  for  the 
slight  broadening  and  displacement  of  bulk-positioned 
atoms  within  the  three  topmost  SiC  bilayers.  Importantly, 
XSW  interpretation  of  the  Si  and  S2  signals  constrains 
these  two  components  into  single,  well-defined  interfacial 
distributions. 

The  XSW-constrained  XRR  result  is  shown  in  Eig.  4(a). 
The  fit  finds  the  EGi,  EG2,  and  EG3  layers  positioned  at 
5.82,  9.17,  and  12.57  A  above  the  terminal  Si  layer  with 
respective  coverages  (©^)  of  0.86,  0.45,  and  0.03  ML.  The 
widths  of  the  graphene  layers  (Table  I),  monotonically 
decrease  as  a  function  of  distance  from  the  interface,  in 
close  agreement  with  trends  observed  in  Refs.  [13,26].  The 
Si  and  S2  layer  positions  converge  to  Zs^  =  2.45  and 

TABLE  1.  XSW-measured  values,  Pj  and  fj,  for  the  jth  chemi¬ 
cal  species.  The  absolute  positions  (zj),  widths  (crj),  and  cover¬ 
ages  (0y)  are  derived  using  combined  x-ray  analysis. 
Uncertainties  on  the  last  significant  figures  are  reported  in 
parentheses. 


j 

Pj 

z/A) 

fj 

(Tj(k) 

^©/MLs) 

fEG3 

12.57 

0.05 

0.03 

EG- 

EG2 

0.42(2) 

9.16 

0.48(5) 

0.11 

0.45 

.EGi 

5.82 

0.21 

0.86 

0.95(5) 

2.4(1) 

0.4(2) 

0.5(1) 

0.70 

52 

0.82(4) 

2.1(1) 

0.9(2) 

0.2(1) 

0.21 

SiBuik 

-0.02(2) 

-0.05(5) 

0.8(1) 

0.3(1) 

0.28 

CBulk 

-0.23(2) 

-0.58(5) 

1.0(1) 

0.1(1) 

0.34 

^Based  on  EG,  1  ML  =  38.2  atoms/nm^. 


215501-3 


PRL  111,  215501  (2013) 


PHYSICAL  REVIEW  LETTERS 


week  ending 
22  NOVEMBER  2013 


Zs^  =  2.13  A,  respectively,  well  within  the  1 -sigma  con¬ 
fidence  levels  derived  from  XSW  analysis.  The  topmost 
bulk  Si  layer  is  —15%  depleted,  which  is  presumably  a 
growth  artifact  from  the  Si  sublimation  and  corroborates 
the  observation  of  near-surface  Si  vacancies  [35].  This  Si 
depletion  limits  the  maximum  number  of  5-^^ -hybridized  C 
atoms  to  —25%  of  the  atoms  within  the  graphenelike 
interfacial  layer.  Finally,  the  C  atoms  in  the  topmost  SiC 
bilayers  are  found  to  relax  outward  by,  on  average,  0.04  A, 
consistent  with  the  XSW  result  (0.05  ±  0.05  A).  With  these 
results  it  is  possible  to  compute  the  XRR-derived  complex 
(0006)  geometrical  structure  factor  for  the  EG  layers. 


3 


I 


(2) 


and  thereby  back-calculate  the  amplitude,  /eg  =  \Seg\  = 
0.49  and  phase.  Peg  =  Arg(5’EG)/27r  =  0.41.  These  XRR 
values  match  closely  with  the  XSW-derived  values  for  EG, 
a  testament  to  the  exactitude  and  self-consistency  of  the 
model  and  analysis. 

The  resultant  chemically  sensitive  atomic  density  profile 
[A^^(z)]  is  shown  in  Eig.  4(b),  and  interfacial  region  detail  is 
shown  in  Eig.  4(c).  Erom  this  construction  we  conclude  that 
the  interfacial  region  contains  a  broad,  C-only  layer  with 
graphenelike  (0.97  ±  0.05  ML)  atomic  density.  Within  this 
layer  there  exist  two  chemically  distinct  C  species.  The  S2 
species,  which  has  the  highest  BE  [Eig.  2(c)],  accounts 
for  —25%  of  the  interfacial  C  and  is  located  2.1  ±0.1  A 
above  the  terminal  Si  layer.  This  displacement  is  only  0.2  A 
larger  than  the  nominal  Si-C  bond  distance  of  1.9  A  in  SiC, 
indicating  the  partial  sp^  hybridization  of  the  interfacial 
layer  with  substrate.  The  narrower  distribution  width  of 
the  S2  species  (cr^^  =  0.2  ±  0. 1  A)  closely  agrees  with  that 
of  the  topmost  Si  layer  (0.3  ±0.1  A),  suggesting  strong 


FIG.  4  (color  online),  (a)  XSW-constrained  x-ray  reflectivity 
results  for  1.3  ML  EG/SiC(0001).  The  best  fit  (black)  to  the  data 
(red)  is  shown  together  with  a  simulated  bulk-truncated  SiC 
curve  (gray),  (b)  The  atomic  density  map  derived  from  combined 
XSW-XPS  and  XRR  analysis,  (c)  The  interfacial  detail  shows 
two  interfacial  components.  Si  and  S2,  separated  by  —0.35  A,  to 
form  a  —0.9  A  wide,  low-areal  density  interface  layer. 


coherency  between  the  two  layers.  Furthermore,  the  cover¬ 
age  ratio  for  the  topmost  Si  layer  to  5*2  is  —1:1  (Table  I). 
These  findings  imply  that  each  Si  atom  in  the  topmost 
SiC  bilayer  is  chemically  bonded  to  an  interfacial  S2  atom 
and  substantiates  the  claim  that  there  exist  essentially  no 
unsaturated  Si  dangling  bonds  at  the  interface  [25]. 

The  Si  species,  which  accounts  for  the  other  75%  of  the 
interfacial  C,  is  positioned  2.4  ±  0. 1  A  above  the  topmost 
Si  layer  and  possesses  a  significantly  larger  distribution 
width  (as^  —  0.5  A)  as  compared  to  the  S2  species.  This 
broadening  reduces  the  maximum  areal  density  at  the  center 
of  the  interfacial  layer  and  is  consistent  with  the  observa¬ 
tions  of  a  highly  corrugated  interfacial  layer  [2,14,22,30]  as 
well  as  the  reduced  effective  interfacial  atomic  density,  as 
seen  with  cross-sectional  transmission  electron  microscopy 
[27,28].  This  broadening  also  points  to  multiple  degrees  of 
strain  in  the  bonding  configurations  for  the  carbons, 
which  is  consistent  with  our  finding  of  a  slightly  broader 
spectral  width  for  the  fitted  Si  XPS  peak  [32]. 

Our  analysis  yields  no  evidence  for  a  number  of  other 
proposed  interfacial  structures.  Specifically,  Refs.  [21,26] 
suggest  the  presence  of  large  populations  of  interfacial  Si 
atoms.  These  models  prove  inconsistent  with  our  combined 
x-ray  analysis.  First,  we  find  no  spectroscopic  evidence  of 
nonbulk  Si  species,  particularly  from  Si  tetramers  or  other 
Si  adatoms,  which  would  exhibit  negatively  core-shifted  Si 
1^  components  with  respect  to  the  SiC  signal  due  to  the 
presence  of  Si-Si  bonds.  Instead,  the  depletion  of  Si  in  the 
topmost  SiC  bilayers  suggests  the  presence  of  negatively 
charged  Si  vacancies  [36],  which  may  play  a  role  in  the 
n-type  doping  of  the  graphene  layers.  With  respect  to  the 
difference  in  interpretation  of  the  XRR  presented  herein 
and  that  of  Ref.  [26],  we  note  that,  because  model-based 
fitting  of  XRR  data  often  leads  to  numerous  statistically 
equivalent  solutions,  unambiguous  analysis  typically 
requires  complimentary  information.  Our  approach  bene¬ 
fits  from  constraints  provided  by  the  structurally  sensitive 
XSW-XPS  measurement  and  finds  that  the  XRR  data  can 
be  well  fit  with  a  C-only  interfacial  layer  with  graphenelike 
density,  therefore  generally  supporting  the  models  of 
Emtsev  [25],  Kim  [14],  and  Varchon  [30].  With  respect 
to  the  discrepancy  in  XPS  peak  fitting  between  this  work 
and  that  of  Ref.  [25],  we  note  the  following:  (i)  the  model 
in  Ref.  [25]  does  not  account  for  the  probable  existence  of 
sub-ML  EG  inclusions  in  their  nominally  zero-layer  gra¬ 
phene,  and  (ii)  there  are  a  large  number  of  potential  influ¬ 
ences  (e.g.,  bond-charge  transfer,  spontaneous  substrate 
polarization  [17],  defect  states  [23],  charge-transfer  doping 
[25],  band  bending  [19],  and  sp^  vs  sp^  hybridization)  on 
the  precise  binding  energies  of  the  interfacial  species, 
making  the  identification  of  chemical  species  from  peak 
shifts  alone  tenuous.  Regardless,  the  binding  energies  that 
we  find  for  5’i(=BEEG  +  0*4eV)  and  5’2(=BEeg  +  LOeV) 
are  in  accordance  with  those  expected  from  sp^-  and 
5-^^ -hybridized  species,  respectively,  and  are  strongly  sup¬ 
ported  by  XSW  analysis. 
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In  summary,  we  have  combined  depth- sensitive  XPS, 
XSW-XPS,  and  XRR  to  create  a  structural  profile  of  the 
EG/SiC(0001)  interface  with  unprecedented  chemical  and 
spatial  resolution.  The  high  degree  of  consistency  between 
the  two  sub-A  resolution  measurements  provides  compel¬ 
ling  evidence  for  the  validity  of  our  structural  model  and 
overall  methodology.  The  interfacial  layer,  which  has  an 
influential  role  in  the  electronic  behavior  of  EG/SiC(0001), 
and,  consequently,  in  graphene-based  devices,  consists 
solely  of  C,  possessing  two  distinct  chemical  species 
located  at  2.1  and  2.4  A  above  the  topmost  Si  layer  of  the 
SiC.  Our  results  support  the  strongly  interacting  interfacial 
layer  model  and  rule  out  any  significant  presence  of  Si 
species  in  adatom  structures  or  within  the  interfacial  layer. 
In  perspective,  XSW-XPS  and  XRR  may  prove  an  effective 
tool  for  the  study  of  the  complex  interfaces  that  exist  in 
emerging  2D  crystals  and  heterostructures,  including  inter¬ 
calated,  doped,  and  functionalized  EG/SiC  nanostructures. 
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X-ray  Standing  Wave 

In  the  dipole  approximation  for  the  photoelectric  effect,  the  normalized  E'^dependent 
photoelectron  yield  is: 

Yj(E^)  =  1  +  R(E^)  +  2p Jr(E^)  fjCos[cf>(E^)  -  27rPj])  ,  (SI) 

where  R  is  the  Bragg  reflectivity,  E^  is  the  incident  photon  energy,  p  is  the  polarization  factor, 
and  (j)  is  the  XSW  phase.  By  fitting  Eq.  SI  to  the  photoelectron  yield  data  from  an  atomic  species 
with  a  specific  chemical  state,  j,  one  can  extract  the  Fourier  amplitude  and  phase  fj  and  Pj  to 
resolve  the  chemically  sensitive  atomic  density  profile,  Nj  (z),  which  is  defined  in  the  main  text 
(Eq.  1). 

Experimental: 

The  0.5  and  1.3  ML  samples  were  grown  from  nominally  on-axis  nitrogen-doped  6if-SiC(0001) 
substrates  graphitized  by  direct  current  flashing  in  UHV.  The  0.5  ML  sample  was  processed  with 
a  550°  C  overnight  degas  followed  by  sequential  flashes  of  1000°  C  for  5  minutes,  1 100°  C  for  5 
minutes  and  1200°  C  for  1  minute.  The  1.3  ML  sample  was  processed  using  the  same  degas  and 


1000°  C  and  1100°  C  anneals,  but  were  treated  with  additional  anneals  at  1200°  C,  1250°  C,  and 


1300°  C  for  2,  2,  and  1  minute,  respeetively. 

EG  synthesis  of  the  1.7  ML  sample  was  earned  out  in  a  eommercial  hot-wall 
Aixtron/Epigress  VP508  chemieal  vapor  deposition  reactor.  Prior  to  graphene  growth,  substrates 
underwent  an  in  situ  H2  etch  at  1520°C  for  30  minutes.  After  etching,  H2  was  purged,  and  the 
subsequent  EG  formation  process  was  conducted  under  a  flowing  Ar  ambient  of  10  standard 
liters  per  minute  at  100  mbar  at  1540°  C  for  30  minutes. 

XSW-XPS  measurements  were  performed  in  UHV  (lx  10"'°  Torr)  at  the  ID32  beam 
line[l]  of  the  European  Synchrotron  Radiation  Facility.  The  incident  photon  energy,  was  set 
near  the  SiC(0006)  back-reflection  condition  (Bragg  angle,  6b  ~88°  and  Ey  -2.450  keV)  using  a 
Si(l  1 1)  double  crystal  monochromator.  The  photon  flux  at  the  sample  surface  was  10^^  photonsA 
within  a  0. 1  X  0.4  mm^  spot  size.  Photoelectrons  were  collected  with  a  SPECS-PHOIBOS  225 
electron  analyzer  positioned  with  analyzer  axis  mounted  parallel  to  the  X-ray  polarization 


FIG.  SI:  Experimental  geometries  for  both  (a)  eonventional  XPS  (a  ~78°)  used  for  survey  seans  and  (b)  highly  surfaee-sensitive 
grazing-emission  XPS  (a  ~2°)  used  for  XSW  measurements.  Tuning  the  emission  angle  to  a  ~  2°  improves  surfaee  sensitivity  by 
effeetively  deereasing  the  sampling  depth  of  the  photoeleetrons  originating  from  deep  within  the  erystal  as  eompared  to  those 
nearer  to  the  surfaee.  The  effeetive  sampling  depth  is  Ag  ~IMFP  Sin(a) 


LEED: 


FIG.  S2:  LEED  patterns  for  both  1.3  ME  UHV-grown  (a)  and  1.7  ME  Ar-grown  (b)  EG/SiC(0001).  Each  image  shows  the 
typical  pattern  with  bright  1x1  EG  (red  arrow)  and  1x1  SiC  (white  arrow)  spots.  The  spots  arranged  in  a  hexagon  about  the 
EG  spots  are  due  to  the  bVSxbVs  i?30°  reconstructed  interfacial  layer.  Ar-growth  resulted  in  larger  surface  domains, 
subsequently  resulting  in  the  sharper  EEED  pattern  in  (b). 

XPS  Survey: 


Figure  S3:  Survey  spectra  for  1.3  ME  UHV-grown  EG/SiC(0001)  (blue)  and  1.7  ME  Ar-grown  EG/SiC(0001)  (red).  Spectra 
were  acquired  with  a  photoemission  angle  a  -78°  (Fig.  SI  (a))  and  using  incident  beam  energies  of  2.450  and  2.465  keV, 
respectively.  The  inset  shows  a  weak  0  1s  signal  present  in  the  Ar-grown  spectrum  associated  with  a  small  amount  of  silicon 
oxide  near-surface  contamination.  Oxide  surface  contamination  is  estimated  in  the  text. 


Binding  Energy  (eV) 


FIG.  S4:  Overlay  of  Si  Is  spectra  from  1.3  ML  UHV-grown  EG/SiC(0001)  taken  at  emission  angles  a  -78°  (blue)  and  a  -2° 
(red).  The  peak  width  broadens  by  -20%  when  measured  using  the  a  =  2°  geometry,  indicating  increased  spectral  contribution 
from  strained  surface  Si  species.  The  difference  trace  is  shown  in  black. 

direction  in  order  to  minimize  the  influence  of  non-dipole  contributions  to  the  photoelectron 
yield  [2,  3].  The  FWHM  total  energy  resolution  of  the  photoelectron  spectra  was  ~0.60  eV, 
which  accounts  for  the  FWHM  incident  beam  bandwidth  of  0.34  eV. 

X-ray  reflectivity  measurements  were  performed  in  ambient  at  the  Advanced  Photon 
Source,  Dupont-Northwestern-Dow  Collaborative  Access  Team  5ID-C  station  using  Ey  =  17.0 
keV  X-rays  collimated  to  a  0. 1x2.0  mm^  spot  size  with  a  flux  of  ~5xl0*'  photonsA.  The  reflected 
intensity  at  the  specular  condition  was  measured  using  an  area  detector  [4,  5].  Below  -0.5  A'\ 
the  finite  surface  domain  size  of  UHV-grown  samples  resulted  in  significant  transverse 
broadening  of  the  specular  rod,  which  inhibited  accurate  integration  of  the  XRR  signal. 

Peak  Fitting: 

For  all  samples,  the  CbuIIc,  Si  and  S2  peaks  are  fit  using  either  pseudo- Voigt  functions  or  with  a 
summation  of  Gaussian  and  Lorentzian  lineshapes  (SGL): 
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where  the  components  are  weighted  by  factor  rj  and  have  common  positions  xo,  and  widths  F.  To 
account  for  the  slight  asymmetry  in  the  Si  Is  peak,  we  used  a  modified  a-SGL  function  [6].  To 
account  for  the  metallic  nature  of  the  EG,  the  peak  is  fit  with  a  Gaussian-broadened  Doniach- 
Sunjic  [7]  profile  (DS): 


Table  SI:  Fitting  parameters  for  C  Is  and  Si  Is  spectra  from  EG/SiC(0001)  samples.  SGL  denotes  a  summation 
Gaussian-Lorentzian,  a-SGL  denotes  an  asymmetric  SGL,  with  asymmetry  factors  a  and  b.  DS  represents  a 
Doniach-Sunjic  curve  with  asymmetry  factor  8. 


0.5  ML  UHV-grown  EG/SiC(0001) 


Si  Is 

C  Is 

Component 

Bulk  Si 

SiOx 

BulkC 

EG 

Si 

S2 

Lineshape 

a-SGL 

SGL 

SGL 

DS 

SGL 

SGL 

e  or  a/b 

0.25,0.09 

- 

- 

0.105 

- 

- 

0.55 

0.25 

0.20 

- 

0.10 

0.20 

Eb 

1841.70 

1844.40 

283.80 

284.80 

285.15 

285.75 

FWHM 

1.25 

2.05 

0.85 

0.70 

1.15 

1.00 

1.3  ML  UHV-grown  EG/SiC(0001) 


Si  Is 

C  Is 

Component 

Bulk  Si 

SiOx 

BulkC 

EG 

Si 

S2 

Lineshape 

a-SGL 

SGL 

SGL 

DS 

SGL 

SGL 

e  or  a/b 

0.25,0.09 

- 

- 

0.105 

- 

- 

0.55 

- 

0.20 

- 

0.10 

0.20 

Eb  (eV) 

1841.65 

- 

283.80 

284.75 

285.10 

285.75 

FWHM  (eV) 

1.12 

- 

0.90 

0.70 

1.1 

1.0 

1.7 

ML  Furnace-grown  EG/SiC(0001) 

Si  Is 

C  Is 

Component 

Bulk  Si 

SiOx 

BulkC 

EG 

Si 

S2 

Lineshape 

a-SGL 

SGL 

SGL 

DS 

SGL 

SGL 

e  or  a/b 

0.25,0.09 

- 

- 

0.105 

- 

- 

0.55 

0 

0.20 

- 

0.10 

0.20 

Eb  (eV) 

1841.65 

1844.10 

283.90 

284.80 

285.10 

285.75 

FWHM  (eV) 

1.1 

1.95 

0.85 

0.64 

1.0 

0.9 

Eq.  S2 


DS{x\z,f,  Xo) 


Cos  ^  +  (1  —  £)Tan  ^ 

[F2  +  (x-Xo)2](l-^)/2 


with  asymmetry  factor  e,  position  xq,  and  width  F. 


All  spectra  were  fit  using  a  Shirley  background  [8].  The  asymmetry  value  8  for  the  EG 
peak  was  set  to  0.105,  consistent  with  observations  from  EG  from  H-intercalated  EG/SiC(0001) 
[9],  Fit  parameters  for  the  spectra  in  Fig.  1  and  S8  are  provided  in  Table  SI. 

XSW  Analysis  Using  Conventional  C  Is  and  Si  Is  Peak-fitting  Models. 

The  peak-fitting  models  used  to  analyze  the  data  in  the  main  text  (and  summarized  in 
Table  SI)  differ  substantially  from  those  typically  employed  [10,  11].  However,  we  note  that  the 
C  Is  data  from  nominally  zero-layer  graphene  presented  in  Ref  [10]  can  be  well  fit  by 
accounting  for  a  small  amount  of  graphene  coverage  and  inverting  the  SiiSi  intensity  ratio  (Fig. 
S5).  The  presence  of  such  graphene  inclusions  on  step  edges  have  been  thus  far  unavoidable 
during  the  production  of  nominally  zero-layer  graphene,  and  are  observed  even  on  the  highest- 
quality  samples  grown  using  state-of-the-art  processes  in  Ar  atmosphere  [12-14].  We  also 
observe  that  the  data  presented  in  the  Ref  [10]  was  acquired  prior  to  the  development  of  more 
well-controlled,  homogenous  EG/SiC(0001)  produced  by  Ar  anneal  [15],  increasing  the 
likelihood  of  relatively  high  EG  coverage  on  the  samples  presented  in  that  work.  It  is  therefore 
likely  that  the  spectra  presented  in  Ref  [10]  should  be  fit  accounting  for  contributions  from  EG 
layers,  as  is  presented  in  Fig.  S5. 
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Binding  Energy 

FIG.  S5:  Fits  to  data  from  nominal  zero-layer  graphene  on  SiC(OOOl)  from  Ref  [10],  with  Ceuik,  EG,  Si,  and  S2  eomponents  in 
blue,  green,  red,  and  brown,  respeetively.  The  speetra  are  fit  aeeounting  for  a  -15%  eoverage  of  EG.  In  eontrast  to  Ref  [10],  the 
81:82  peak  intensity  ratio  is  essentially  inverted. 
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FIG.  86:  XP8  peak- fitting  and  subsequent  X8W  data  using  the  fit  parameters  suggested  by  Emtsev  et  al.  and  Riedl.  The  data 
presented  in  these  figures  are  the  same  as  that  presented  in  Figs.  2(e)  and  (d).  (a)  The  C  Is  data  were  fit  with  three  peaks  for  X8W 
analysis  beeause  the  81  and  EG  peaks  were  statistieally  inseparable,  (b)  The  8i  Is  peak  fitting  model  shown  here  aeeounts  for  the 
possible  presenee  of  distinet  6R3  and  defeet  related  eore-shifted  eomponents.  In  (e)  and  (d)  the  X8W  results  eorresponding  to  the 
peak  fitting  models  in  (a)  and  (b),  respeetively.  XP8  yield  eurves  are  offset  on  the  y-axis  for  elarity. 


Figure  S7:  Goodness-of-fit  maps  for  the  Si-Si  distance  with  fixed  S2  position  using  traditional  [10]  peak-fitting 
models.  There  exist  two  local  minima,  indicating  possible  solutions  at  Si-Si  =  0.9  A,  and  Si-Si  =  2.4  A,  but  both 
solutions  lack  realistic  physical  interpretation. 


Similarly,  Riedl  et  al.  propose  that  the  Si  2p  spectrum  be  fit  with  a  3 -peak  model  due  to 
the  presence  of  6R3  and  “defecf’  species  [11].  We  note,  however,  that  even  a  moderate 
population  of  surface- specific  species  would  presumably  dominate  the  spectrum  in  a  fashion 
similar  to  that  observed  for  the  C  Is  spectrum,  while  the  observed  increase  in  the  relative 
intensity  on  the  wings  of  the  Si  Is  peak  are  marginal  (~x3)  when  the  measurement  taken  in  the  a 
~2°  geometry  as  compared  to  that  taken  in  the  conventional  geometry  (Fig.  S4).  Therefore,  we 
advocate  that  the  increased  intensity  at  the  wings  does  not  arise  from  distinct  Si  species  with 
discrete,  core-shifted  components,  but  is  rather  caused  by  the  distribution  of  bond  angles  and 
bond  lengths  present  at  the  strained  interface.  Finally,  we  note  that  in  the  work  of  Riedl  et 
n/.[ll],  both  the  6R3  and  “defect”  components  disappear  upon  H-decoupling  of  the  buffer  layer 
from  the  substrate,  which  is  consistent  with  our  interpretation  of  spectral  broadening  due  to 


strained  Si  species. 


Table  S2:  XSW  results  based  on  XPS  peak- fitting  models  from  Refs.  [10]  and  [11].  Reported  uneertainties  are  1 -sigma 
eonfidenee  bounds. 


Z/(A) 

fi 

o;(A) 

BulkC 

1.30 

0.76+0.03 

2.39+0.13 

0.85+0.1 

0.23^H^ 

C/3 

EG  +  Si 

1.53 

0.39+0.03 

N/A 

0.22+0.03 

N/A 

u 

S2 

2.94 

0.83+0.04 

2.09F0.10 

0.65+0.15 

n  Qy+o.io 

Bulk  Si 

4.85 

1.00+0.02 

2.52+0.05 

0.88+0.08 

0.21^H^ 

c/3 

Si6«5 

2.15 

0.95+0.06 

2.42F0.13 

0.9F0.3 

xn 

Sidef 

0.79 

0.97+0.11 

2.44F0.25 

0.6+0.3 

Regardless,  we  are  not  limited  in  the  analysis  of  a  single  peak  fitting  model,  and  therefore 
provide  alternate  XSW  analysis  using  the  XPS  models  offered  by  Emtsev  and  Riedl.  The  C  Is 
and  Si  Is  spectra  fit  according  to  these  models  in  shown  in  Fig.  S6.  Both  spectra  are  well  fit 
based  on  the  literature  values.  We  note  that  in  Emtsev’s  model,  the  Si  and  EG  components  differ 
in  binding  energy  by  only  ~0.1  eV.  This  results  in  a  high  degree  of  covariance  of  peak  fitting 
parameters  for  these  two  species,  greatly  complicating  the  XSW  analysis.  The  S2  component,  on 
the  other  hand,  is  practically  isolated  (core-shifted  by  +1.1  eV)  and  can  be  analyzed  in  a 
straightforward  manner.  The  XSW  data  and  fits  are  presented  in  Figs.  S6(b)  and  (c),  and  the 
results  are  summarized  in  Table  S2.  Because  direct  analysis  of  the  EG  +  Si  component  proves 
impractical  (the  one  measured  Fourier  component  possesses  as  many  as  4  distinct  contributions), 
we  explore  Emtsev’s  model  by  constraining  the  S2  species  within  its  XSW-derived  1-sigma 
confidence  window  (2.0  A  <  zs2  <  2.2  A,  see  Table  S2)  and  mapping  the  goodness-of-fit  x^and 

R-factor  as  a  function  of  zsjbetween  0.5  >  zsj  <  3.0.  The  resulting  x^and  R-factor  maps  are  shown 

in  Figure  S7.  There  exist  two  distinct  local  minima  in  the  map,  indicating  potential  solutions  for 
Emtsev’s  model  at  zsj  -0.9  A  zsj  -2.4  A.  The  zsj  -0.9  A  solution  would  indicate  Si-C  and 


graphene-like  C-C  bonding  distances  that  are  incompatible  with  the  interpretation  of  a  partially- 


bound  graphene-like  interfacial  layer.  The  zs^  ~2.4  A  solution  produces  a  structure  largely 

similar  to  the  one  presented  in  the  main  text,  but  places  the  Si-C-Cs  bonded  atoms  at  distances 
much  larger  than  typical  Si-C  bond  lengths.  Furthermore,  this  model  places  the  Si-C-Cs  bonded 
C  atoms  in  a  highly  unphysical  bonding  geometry,  at  0.3  A  above  the  atoms  in  a  graphene-like 
configuration. 

The  XSW  modulations  resulting  from  the  fitting  of  the  Si  Is  spectra  with  distinct,  core- 
shifted,  non-bulk-like  components  are  shown  in  Figure  S6(d).  The  Fourier  amplitudes  and  phases 
for  these  three  components  are  practically  indistinguishable  within  error  (Table  S2),  indicating 
all  species  have  similar  positions  and  distributions  with  respect  to  the  substrate  lattice.  While  it  is 
possible  that  this  result  indicates  the  existence  of  small  populations  (<5%  at  surface)  of  distinct 
chemical  species  positioned  in  bulk-like  SiC  positions,  we  would  expect  some  non-bulk  like 
XSW  modulation  of  the  Si  Is  spectral  components  if  any  Si-Si  bonds  were  present,  and  none  is 
observed.  An  alternate  explanation,  in  which  interfacial  strain  causes  increased  variance  in  the 
distribution  of  Si-C  bond  lengths  and  angles  of  the  topmost  Si  layers,  thereby  broadening  the  Si 
Is  spectral  distribution,  is  more  plausible. 

Analysis  of  0.5  ML  UHV-grown  1.7  ML  Ar-grown  EG/SiC(0001): 

C  Is  and  Si  Is  spectra  taken  with  emission  angle  a  =  2°  are  shown  in  Fig.  S8  for  both  0.5  ML 
UHV-grown  [S8(a)-(b)]  and  1.7  ML  Ar-grown  [S8(c)-(d)]  EG/SiC(0001).  Both  samples  exhibit 
C  Is  spectra  typical  of  EG/SiC(0001),  but  both  Si  Is  spectra  also  possess  a  strong  high-BE 
component  consistent  with  SiOx  chemical  species  [16].  This  signal  is  only  discernible  when 
using  the  a  =  2°  geometry,  indicating  that  it  is  associated  with  a  surface  oxide  species.  Due  to  the 
relatively  low  spectral  resolution  of  our  measurement,  we  do  not  attempt  to  identify  sub-oxide 
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Figure  S8:  C  Is  and  Si  Is  spectra  from  0.5  ML  (a)-(b)  UHV-grown  EG/SiC(0001)  and  (c)-(d)  1.7  ML  Ar-grown  EG/SiC(0001). 


Table  S3:  XSW  results  for  0.5  UHV-grown  and  1.7  ML  Ar-grown  graphene. 


UHV-grown  0.5  ML  EG/SiC(0001) 


Component,  s 

x' 

Ps 

/s 

a*  (A) 

Bulk  Si 

2.12 

0.9910.01 

2.4910.03 

0.9010.03 

SiOx 

2.03 

N/A 

N/A 

0.110.1 

N/A 

BulkC 

6.22 

0.7410.02 

1.8710.05 

0.9510.08 

n  -1  +0.07 
^■-*--0.1 

Si 

1.36 

0.910.1 

2.310.2 

0.310.2 

o.6lg:i 

S2 

1.21 

0.810.05 

2.010.1 

0.910.2 

o.i8^«ii 

EG 

1.45 

0.2710.04 

N/A 

0.7010.3 

N/A 

Ar-grown  1.7  ML  EG/SiC(0001) 


Component,  s 

x' 

Ps 

2s(A) 

fs 

a*  (A) 

Bulk  Si 

0.33 

1.0010.01 

2.5210.03 

0.9410.08 

'X  +0.07 
'-'■-*--0.1 

SiOx 

1.09 

N/A 

N/A 

n  n+0.2 

N/A 

BulkC 

3.37 

0.7410.02 

1.8710.05 

0.8610.07 

n  99  +  0.06 
U.ZZ.-0.06 

Si 

3.11 

1.0010.05 

2.5210.13 

0.310.1 

o.6^g:J5 

S2 

0.86 

0.8410.03 

2.1210.07 

1.0010.14 

rt  'i  +0.07 
'-'■-*--0.1 

EG 

1.21 

0.3410.04 

N/A 

0.5210.09 

N/A 

signals  within  the  Si  Is,  as  have  been  observed  previously  [16].  We  do  note,  however,  that  in  the 
case  of  the  1.7  ML  sample,  which  has  a  weaker  oxide  signal,  the  SiOx  core-level  shift  is  0.3  eV 
less  than  that  of  the  0.5  ML  sample,  suggesting  increased  contribution  from  lower-BE  sub-oxide 
species.  The  oxygen  concentration  can  be  estimated  by  comparing  the  integrated  O  Is  signal 
intensity  to  that  of  the  EG  C  Is  signal.  The  density  of  C  in  EG  form  is  estimated  from  the  unit 
cell  definition  to  be  38.2  C/nm  ,  and  therefore,  by  correcting  for  the  relative  C  and  O 
photoionization  cross-sections  at  ~2.5  keV,  we  estimate  a  concentration  of  ~2  0/nm  for  the  1.7 
ML  Ar-grown  sample,  and  ~6  0/nm  for  the  0.5  ML  UHV-grown  sample.  For  the  Ar-grown 
sample,  the  oxide  formation  may  be  caused  by  the  presence  of  trace  O2  or  H2O  in  the  chamber 
during  growth  or  cool  down.  In  the  case  of  the  UHV-grown  sample  the  oxide  formation  may 
arise  due  to  exposure  of  the  more  reactive  sub-ML  sample,  which  has  large  regions  of  exposed 
zero-layer  graphene,  to  air. 

XSW  results  for  both  samples  are  shown  in  Fig.  S9,  and  the  results  are  summarized  in 
Table  S3.  The  most  notable  result  is  that  for  the  Sisuik,  Csuik,  Si  and  S2  values  are  essentially 
identical  between  all  measured  samples.  Considering  that  these  samples  were  made  using 
different  growth  methods,  in  different  laboratories,  and  possess  differing  EG  and  oxide  coverage, 
these  nearly  identical  results  are  a  testament  to  the  ubiquity  of  the  interfacial  structure.  The  SiOx 
components  are  randomly  distributed,  indicating  that  there  exists  a  thick  or  broadly  distributed 
(as  opposed  to  single  oxide  monolayer)  region  of  silicon  oxide  near  the  SiC  surface.  The  0.5  ML 
UHV-grown  EG  result  is  interesting  as  it  approaches  the  value  expected  for  mainly  monolayer 
graphene  (at  zegj  ~5.8  A,  an  isolated  monolayer  would  give  a  coherent  position  of  Pegj  =  0.30). 

This,  along  with  the  relatively  high  measured  coherent  fraction  if^o  =  0.7+0. 3),  indicates  that 
most  of  the  EG  signal  originates  from  monolayer  graphene.  Contribution  from  the  second  layer 


of  EG,  EG2  (at  ZEG2  would,  in  principle,  shift  the  Fourier  phase  positive  and  reduce  the 

coherent  fraction,  as  is  observed  for  samples  with  higher  EG  coverage.  It  should  be  noted  that  the 
XSW-XPS  data  for  the  0.5  ML  sample  possess  poorer  counting  statistics,  which  is  reflected  in 
the  generally  larger  degree  of  scatter  in  the  data  and  uncertainty  values  in  the  extracted  XSW 
results. 


The  1.7  ML  Ar-grown  sample  serves  to  validate  the  XSW-XRR  analysis  performed  in  the 
main  text.  Qualitatively,  the  XSW  values  for  Si  and  S2  agree  extremely  well  with  those  in  the 
main  text  (within  1 -sigma  in  both  cases),  indicating  that  the  two  samples  have  essentially 
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FIG.  S7:  XSW  results  for  (a)  C  Is  0.5  ML  UHV-grown.  (b)  Si  Is  0.5  ML  UHV-grown.  (c)  C  Is  1.7  ML  furnace-grown,  (d)  Si  Is 
1.7  ML  UHV-grown  EG/SiC(0001)  samples. 


FIG  S8:  XRR  analysis  for  the  1.7ML  Furnace-grown  EG/SiC(0001)  sample.  Results  are  highly  comparable  to  those  found  in  the 
main  text,  apart  from  a  higher  degree  of  graphene  coverage. 

identical  interfacial  structure.  The  differing  XSW  result  for  the  EG  component  is  due  to  the 
varying  amount  of  EG  coverage  on  the  samples.  The  XRR  data  and  XSW-constrained  best-fit 
result  is  shown  in  Fig.  SIO.  We  find  that  the  values  zegj  ~5.80  A,  zeg2  ~9.15  A,  and  zegj  -12.55 

A  yield  y^=  9.69  and  f?-factor  =  0.085,  comparable  to  values  found  for  1.3  ML  EG/SiC(0001)  of 
%  =7.19  and  R-factor  =  0.079  .  The  slightly  poorer  quality  of  the  fit  for  the  Ar-grown  sample 
may  be  attributed  to  the  inability  to  accurately  model  the  contribution  of  the  SiOx  in  the  electron 
density  profile.  Relative  layer  coverages  were  cegi  -LOO  ML,  ceg2  -0.55  ML,  and  cegj  -0.10 

ML,  summing  to  a  total  coverage  of  1.7  ML.  The  back-calculation  of  the  XSW  values  gives ^eg  = 
0.38  and  Peg  =  0.44,  near  the  1-sigma  of  the  uncertainty  limits  for  XSW  results  reported  in  Table 


S3.  In  all,  the  consistency  across  sample  coverage,  production  method,  and  sample  morphology 
provides  strong  support  for  conclusions  stated  in  the  main  text. 
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